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We present a perturbative approach within the scope of 
Kohn-Sham density functional theory (DFT). The method is 
based on the exact exchange-only optimized effective poten- 
tial method, and correlation is included via perturbation ex- 
pansion using Rayleigh-Schrodinger (RS) perturbation theory 
(PT). The correlation potential is constructed from M0ller- 
Plesset formulation of RS calculations. This naturally leads 
to a new iterative scheme when finite order perturbation the- 
ory is employed. The new iterative procedure can be taken as 
a self-consistent parameter-free DFT PT, and, as such, pro- 
vides correlation energy which is explicitly functional of the 
self-consistent orbitals, eigenvalues and also of the exchange- 
correlation potential. As a demonstration, terms up to the 
second order in the PT expansion are considered and the 
explicit formula for the second order correlation potential is 
given as well. 



I. INTRODUCTION 

Kohn-Sham (KS) density functional theory (DFT) is a 
formally exact treatment of the many-body problem [Q . 
Kohn and his co-workers showed that the exact equations 
of quantum mechanics, taking into account all the elec- 
trons in a molecule (or solid) can be replaced by simpler 
equations in which only the density of electrons at each 
point in space enters Remarkable success has been 
achieved in the last decades in finding the " missing link" 
of KS-DFT, the accurate exchange-correlation energy 
functional E xc [p]. The possibility of further improve- 
ments of model functionals, however, seems to be rather 
limited within the scope of the so-called second gener- 
ation of DFT, in which the kinetic energy is expressed 
in terms of orbitals, while the exchange-correlation en- 
ergy E xc as an electron density functional. Furthermore, 
"state-of-the art" exchange-correlation energy function- 
als E xc [p] and potentials v xc (r), such as the generalized 
gradient approximation (GGA) Q suffer from serious 
problems since they do not exhibit the correct — £ long- 
range and quadratic behaviour at small r HQ. It be- 
came also apparent that GGA is not able to reproduce 
exchange-correlation energies and potentials simultane- 
ously with the accuracy provided by the energy function- 
als Q. All this indicates that although density gradient 
treatments have had some success, they are not sufH- 
ciantly well developed at present to provide a final answer 
to the problem of determining the exchange-correlation 



contribution to the ground-state energy functional or to 
the one-body potential ||. 

A new generation of DFT, the optimized effective po- 
tential method (OEP) opens a new horizon towards 
fully local and exact mean-field theory, where not only 
the kinetic energy but also the exchange-correlation po- 
tential and the exchange-correlation energy are treated 
at the orbital dependent level of theory. The efficient 
and highly accurate KLI approximation to OEP pro- 
vides an efficient computational tool for the scientific 
community. OEP and KLI build up correct — £ assymp- 
totics and the highest occupied orbital energy satisfies 
Koopmans-theorem and also provides self-interaction free 
exchange-correlation potential v xc (r) S. However, the 
theory is exact at the exchange-only level and electron 
correlation can only be added via approximate corre- 
lation energy functionals |]]. It must be emphasized 
that the correlation energy functionals neither in their 
local nor in gradient corrected form are likely to work in 
combination with exact exchange. This is because ba- 
sically none of the common correlation functionals have 
a long-range component in the corresponding correlation 
hole (the combined exchange-correlation hole is typically 
short-ranged). The Colle-Salvetti gradient corrected cor- 
relation functional |^], which provided excellent results 
for atoms, performed rather badly for molecules due to 
the abovementioned reason H . Approximate correlation 
functionals, which are derived from the homogeneous 
or inhomogeneous electron gas model or obtained non- 
empirically from sum-rule conditions, exhibit incorrect 
long-range tail (only dynamical correlation is accounted 
for). Those functionals also build up improper local 
behaviour in atoms and molecules and their success is 
mainly due to the cancellation of errors with differrent 
sign fp|JlC|] . Therefore, the correlation effect can only be 
calculated in a more complicated way than in the second 
generation of DFT in those methods which work with 
exact exchange (e.g. HF and OEP) pfj| . 

A few extensions of the OEP method have been made, 
which allow the treatment of electron correlation in a lim- 
ited way. In recent years much effort has been devoted to 
the formulation of perturbation theory (PT) on KS basis. 
Recently Rayleigh-Schrodinger PT has been employed 
for the DFT correlation energy using explicit coupling- 
constant dependence in the model Hamiltonian within 
the KS density functional picture pi . The Gorling- 
Levy scheme and OEP provide the exact treatment of 
the exchange-only KS problem and also provide an alter- 
native to the Hartree-Fock theory. Casida, on the other 
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hand, employed an approximate perturbative expression 
for the ground state to include correlation in OEP us- 
ing a formalism based on the work of Sham and Schlutcr 
In a recent review of OEP, Grabo et al. proposed 
a scheme which also makes use of many-body pertur- 
bation theory (MBPT) on exchange-only OEP reference 
state using Green's functions and the Dyson equation 
fl4|| . Very recently, Engel et al. worked out the rela- 
tivistic generalization of OEP and used M0ller-Plesset 
based correlation |l5|]. Aashamar et al. gave the exten- 
sion of OEP approach which, permits the treatment of 
correlation in a limited way at the multiconfigurational 
self-consistent field (MCSCF) level of theory @. How- 
ever, this approach met little practicle relevance being 
extremely time-consuming and also difficult to relate to 
the single determinent nature of KS DFT. 

The density functional perturbation theory worked out 
by X. Gonze et al. which can be taken as a variation- 
perturbation procedure. In this as well as in other the- 
ories the self-consistent procedure led to the Kohn- 
Sham wave function being dependent on the perturbation 
theory. Thus the perturbation comes into play at an ear- 
lier stage than in conventional quantum chemical pertur- 
bation theory (e.g. M0ller-Plesset perturbation theory). 
Also, Holas and March expressed the exact exchange- 
correlation potential in terms of first- and second or- 
der density matrices using the Gorling-Levy perturbation 
theory |nj . All these methods appear somewhat complex 
in detail and it is presently difficult to predict which one 
will prove the most convenient for practical application 
and implementation fTsf] . 

The Kohn-Sham formalism for ground states is based 
on the noninteracting Schrodinger equation 



\T- 



(i) 



the KS equation, where v s is the A^-electron opera- 
tor which corresponds to a local multiplicative potential 
v s (r) as a consequence of the Hohenberg-Kohn theorem 
The one-body KS potential v s (r) is determined, up to 
an additive constant, by the requirement that the ground 
state of the KS Hamiltonian operator T+ii s , the KS wave 
function <&, yields the same electron density p(r) as the 
ground state of the corresponding interacting real system. 
$ can be composed of N spin orbitals Ui as a single Slater 
determinant for nondegenerate systems. Therefore, the 
following spinrestricted Kohn-Sham single-particle ap- 
proach is considered for fermionic systems (we have omit- 
ted spin-dependency, but the extension of the theory is 
straightforward for that case, throughout a.u. is used): 



[-^V 2 + v s (y)]u 1 (y) = ejUi(r) 



(2) 



where Ui(r) and ti are the single particle orbitals and 
eigenvalues for a fermionic system and v s (r) is the effec- 
tive Kohn-Sham single particle potential 



where v ex t is the potential external to the electronic sys- 
tem that includes the one created by nuclei, vh is the 
Coulomb potential due to the classical electron-electron 
repulsion and nonclassical term, the exchange- 

correlation contribution to the the KS one-body poten- 
tial. 

In this article we would like to study a more general 
class of KS potentials. In a general sense, the exchange- 
correlation part of the KS potential is not a pure density 
functional but rather a complicated functional of eigen- 
values and of single particle orbitals [M, 



Vs([p],r) =v s ([p,{ui},{ei}],r) 
= v ext ([p],r) +v H {[p], r ) 
+v xc ([p,{ui},{ei}},r) 



(4) 



In addition to this, in the section V we will show that 
the exchange-correlation energy E xc is even more com- 
plicated energy functional being the explicit functional 
of the exchange-correlation potential v xc (r) as well, 



E r 



(5) 



The same conclusion is derived by others |14| . As demon- 
strated by Gorling and Levy [^T) , it is not necesseraly im- 
portant to know how the exchange-correlation potential, 
based on the eigensolutions of the KS equations, depends 
on the density. Those potentials are always implicitly 
functionals of the density. 

To derive v xc formally, we write E xc rigorously 



Exjp] = 



1 



dX / dr'dr 



P (r)p(r')[ff A ([p];r,r')-l] 



(G) 



where g x ([p];r,r') can be interpreted as the pair corre- 
lation function of the fictitious system with interaction 
strength parameter (coupling constant) < A < 1 and a 
ground state density which is independent of A (T||4|] . The 
exchange-correlation potential v xc is formally defined as 
the functional derivative 



(r,[p}) = 



5E xc [p] 
Sp(r) ■ 



(7) 



Because of the Hohenberg-Kohn theorem, there exists a 
one-to-one mapping between the single particle KS po- 
tentials v s (r) and the densities p(r) which guarantees 
that the functional derivative given by Eq. (Q) is defined 
Q. However, the functional derivation can not be ob- 
tained explicitly, since one has to derive the following 
expression 



,(r) = v ext (r) + v H (r) + v xc (r), 



(3) 



Jo I J I 1 - 1 I 

I f ^,J' P(r')p(r")6g x ([p];r\r") 
2j drdr |r'-r"| Sp(r) 
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The first term of the potential is the so-called potential 
part of the exchange-correlation energy |Q] , which is iden- 
tical with Slater's potential in the exchange-only (x-only) 
case, i.e. for o A approximated by g A=0 , and which is ex- 
actly known [|23| as an orbital dependent quantity. The 
response like second term contains the fundamental func- 
tional derivative of the pair correlation function which is, 
however, an unknown functional. The functional deriva- 
tive Eq. (0) can be calulated only for approximate func- 
tionals analitically where explicit dependence on the den- 
sity p is known and therefore to be applied only for the 
so-called second generation of DFT jl4|] . In the third gen- 
eration of DFT one uses -E x [{iti}] rather than E x [p] so 
that not only the kinetic energy but the exchange energy 
is expressed as orbital dependent energy functionals. The 
correlation energy must be determined, though as a den- 
sity functional p|,|l4||. The central equation in the third 
generation of DFT is still the KS Eq. (§) At the 

best of our knowledge attempts have not yet been made 
to express the exchange-correlation energy as a unique 
orbital-dependent quantity -EzcK^j}] in OEP. Further- 
more, a general theory to be set up, where both the 
exchange-correlation potential and the energy are com- 
plicated functionals of the self-consistent solutions of the 
eigenvalue problem Eq. (Q) in the spirit of Eqs. (3) and 
(4). 

In this work we give an explicit formulation of a pertur- 
bation theory on top of exchange-only OEP. We derived 
an explicit expression for the correlation potential given 
by Krieger et al. ||. In section II, we give the short 
summary of OEP formalism, which we use extensively in 
the further sections of this article. We give in sections 
III-IV, the potential and response part of v xc (see Eq. 
(7)) exactly. We formulate a self-consistent perturbation 
theory on top of the exchange-only OEP reference state 
(section V) . We would also like to compare our OEP-PT 
scheme with other DFT schemes obtained by perturba- 
tion theory |11 13] in order to point out the similarities 



and differences between them. 



II. OPTIMIZED EFFECTIVE POTENTIAL 
METHOD 

The starting point of the OEP method is the total 
energy functional 

Egf P [{M] =E/ dvu*{v) ( -iv 2 ) Ui (r) 

+ J drp(r)v ext (r) + \ J dvdv' ^£^- + E g EP [{«,}] , 

where, in contrast with ordinary DFT, the exchange- 
correlation energy is an explicit functional of orbitals, 
and, therefore, only implicit functional of the density via 
Eqs. (2)-(3) P,H- i is a collective index for all orbital 



quantum numbers. The local single-particle potential ap- 
pearing in Eqs. (§)-(3) must be the optimized one yield- 
ing orbitals which minimize Ef ) f p [{u i }] so that 



SE°? p [{u t }} 
Sv s (r) 



= 0. 



(10) 



As first pointed out by Perdew and co-workers, Eq. ( pT| ) 
is equivalent to the Hohenberg-Kohn variational principle 

The optimized effective potential method (OEP) is 
given by Talman and Shadwick Jt|] for getting the ex- 
act exchange potential v x (r). Formally, we make use of 
chain rule for the functional derivative Eq. (6) to obtain 



SE KC [{ut\p]}] 



Sp(r) 



(11) 



E 

i=l 



dr , 6E° EP [{u i }]6u i (r') | cc 
Sui(r') 6p(r) 



where are the orbitals which are, however, implic- 

itly functionals of the density. Applying the functional 
chain rule again and after some algebra one can get the 
following integral equation [Bill 24 1 : 



£>(r) J dr / [«° BP (r / )-« i (r / )]G i (r,r / K(r') + c. C . = 

i 

(12) 



where 



Vi(r) 



1 SE xc [{ Ui }} 



<(r) S Ui (v) ' 
and Gi(r,r') is the Greens function 

G t r,r = > - • 

3=£* 



(13) 



(14) 



The integral Eq. (12) is the fundamental expression for 
v xc in OEP instead of Eq. however, there is no known 
analytic form of u xc [{Mi}]. Therefore, only numerical so- 
lutions are available for spherical atoms [|3].]7|, pT| , pii|Jl7| ] and 
for solids [p5|-f27| . These numerical solutions of Eq. jl^ ) 
are confined to the exchange-only OEP and correlation 
has been taken into account only via approximate local 
functionals ^5|. At the best of our knowledge, the only 
possibility to account for the correlation energy exactly 
is the constrained search formulation of KS theory, when 
one makes use of fully correlated density as a reference 
density obtained from e.g. full CI calculations P,|l0|]. Al- 
(9) though such a calculation provides exact v xc , it repre- 
sents an enormous computational task which can hardly 
be carried out for systems of practical interest. Gorling 
and Levy have pointed out that Eq. (|lj) is the special 
case of a more general equation |2s| . 

In the work of Krieger and co-workers, the OEP inte- 
gral equation is analyzed and a simple approximation is 
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made reducing the complexity of the original OEP equa- 
tion significantly, and, at the same time, keeping many of 
the essential properties of OEP unchanged S. Krieger, 
Li and Iafrate gave an exact expression by transforming 
the OEP integral Eq. (7) into a manageable form. They 
have got the following, still exact expression for v xc (r): 



-HF 



J xcj 



drpj(r)vf c (r) 



dv Pj (r)v C j{v). 



occ 



,-,..-'(r) = <4(r) + > J PiW {t OEP_ Vi) (15) 



^ P(r) 

' V[p l (r)Vw,(r)] 



P(r) 



where ff c (r) is the Slater's potential given as the first 
term in Eq. (7) in exchange-only case. However, Slater's 
potential can be generalized to include the potential part 
of the correlation energy density as well The ex- 

change component of uf c 



OCC p 

^ (r) = ~2^^" l(r) ^* (r) / dt 

i,3 



,<(r')^(r') 



r — r 



(16) 



The summation runs over the orbital index for all the 
occupied orbitals up to the highest occupied mth orbital 
(Fermi level) . The function pi (r) is defined by 



Pi(r) 



1 



«i(r) 



dr'[^ p (r')-^(r')]G,(r,r')^(r'), 



(17) 



with the partial density pi(r) — u*(r)ui(r). In practical 
applications the last term in Eq. (15) turns out to be 
quite small in atomic systems and has a small effect only 
on the atomic shell boundaries ||. The neglect of this 
term then leads to the KLI-approximation which has the 
following form after some algebra H. 



v. 



m— 1 / \ m — 1 



P(r) 



3=1 



(18) 



m — 1 



( A X )y ( V xcj ~ V j) = '''r,, 



—KLI — 



3=1 



-HF 



4^ 

i=l 



dvdv 



,Ui(r)^(r')u|(r')^(r) 



(19) 



where m is the Fermi level (the highest occupied one- 
electron energy level) and note that the mth level is ex- 
cluded because ii™ = v^ F This might appear a 
rather crude approximation, but it can be interpreted as 
a mean-field approximation since the neglected terms av- 
eraged over ground-state density vanish. It can be shown 
that the KLI approximation can be derived alternatively 



III. PERTURBATION EXPANSION OF THE 
CORRELATION ENERGY ON TOP OF 
EXCHANGE-ONLY OEP 

Let us turn now to the discussion of how to include 
electron correlation on top of exact exchange-only KS 
theory. Our aim is to solve the exact KS scheme using 
perturbation theory on top of a known first order problem 
(e.g. the x-only OEP). Here consider the Hamiltonian for 
an TV-electron system 



H 



Vee + y^^Vextjn), 



(20) 



where T = X)i=i — 2^1 > V ee = V ee (r) stands for the 
electron repulsion as a local operator and v ext (ri) is the 
external (nuclear) potential of the nuclear frame. 

We consider the problem of improving the exchange- 
only OEP energy of an TV-electron system by means 
of Rayleigh-Schrodinger (RS) perturbation theory. The 
partition of the Hamiltonian is defined H as 



H = Hn 



w 



(21) 



In particular, we treat W as a perturbation to the 
exchange-only KS Hamiltonian H (H >> W). We 
wish to solve the eigenvalue problem 



M„= /^^*,< > i = l,...,m-l 



7=1 



H$ = (H + W)$ = E<£, 



(22) 



which leads to a non-degenerate ground state for a system 
of TV-fermions. $ is the ground state wavefunction and 
E is the expectation value of the Hamiltonian given by 
Eq. ( jgCj ) and (|2l|). Eq. (p2|) describes the corresponding 
KS system, a model system of hypothetical noninteract- 
ing electrons with the same ground-state electron density 
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as the real system [Jl|j5|Jl 1[ ■ However, we only know the 
solution of the unperturbed problem 



H $ s = E Q <S>s 



(23) 



where Ho, $ s are the unperturbed reference Hamiltonian 
and single determinent wavefunction. We assume that 
Eq. (|||) provides a complete set of eigenfunctions {ui} 
with corresponding eigenvalues {e^} as well as the nonde- 
generate ground energy state Eq . In this article we treat 
Ho as the exact exchange-only KS Hamiltonian, which is 
known as the exchange-only OEP Hamiltonian. Actually, 
the solutions of exact x-only OEP are equivalent to the 
solutions of the exact x-only KS equations jlj} . The cen- 
tral idea is now to determine the perturbation operator 
W as a functional derivative with respect to the electron 
density p (Eq. (§). 



W = W(r) 



<y(£7-(4.|H |$.» 



Vee(r) 



5p 



(24) 



Therefore, the perturbation operator W is treated as a 
local potential W(r). The calculated operator W can 
be then substituted back to the Hamiltonian given by 
Eq. ( pl[ ) and the eigenvalue problem (Eq. (p2|)) can be 
solved. The unknown ground state energy E in Eq. (24). 
is expressed by RS perturbation theory. In section V 
we give the details of this iterative scheme, in particu- 
lar when the reference state is the exchange-only OEP 
method. 

We treat v s ( [p] ; r) in this article as an explicit func- 
tional of the eigensolutions {ui, e,} of the single particle 
Hamiltonian (to be interpreted as the KS one). There- 
fore, the unknown potential v xc is implicit functional of 
the ground state density according to the Hohenberg- 
Kohn theorem [Q. Its uniqueness is also guaranteed by 
the Hohenberg-Kohn theorem. v s ([/?];r) in Eq. ph is a 



one-body operator and Eq. (|23j) reduces to one-particle 
equations for KS single particle orbitals Ui (r) as given by 
the noninteracting one-particle Schrodinger Eq. (||) with 
the exact x-only KS one-body operator v^ OEP . The KS 
single determinant, 5> s , in Eq. ( |23| ) is then formed from 
all occupied one-particle KS orbitals, {itj}, i.e., the TV en- 
ergetically lowest solutions to Eq. ( p3|) with N being the 
number of electrons in the system of interest. Eq. (|22| ) re- 
duces to exact one-particle KS equations given by Eq. (^) 
with the exact KS one-body operator v s (r). Further- 
more, the one-body potentials, which correspond to the 
x-only reference state (Vg OEP ) as well as to the exact 
(correlated) case (v^ s ) together with the Hamiltonians 
Ho and H, are unique functionals of the u-representable 
density. However, it must be emphasized that the x-only 
OEP will result in ground state density different from the 
exact KS scheme given in Eq. (^2|). 

The functional derivation given by Eq. (24) can not be 
carried out directly due to the difficulties given in Eqs. 
(7)-(8). The derivative is directly not accessible since E xc 



is known only in terms of the one-particle KS states {ui} 
and the eigenvalues {e^}, and the explicit functional de- 
pendence of the KS eigensolutions on the electron density 
is unknown. OEP provides an alternative way of getting 
the local operator W(r). The solution of Eq. ( p3|) is 
identical with the solution of Eq. (|l^) for exchange-only 
OEP. Therefore, we chose the following partition of H 
according to Eq. (|22|), 



(hi+vrn 



and 



W(r) = V ee (r) 



(25) 



(26) 



where hi is the one-electronic Hamiltonian which con- 
tains the kinetic operator and the operator of the exter- 
nal (nuclear) potential (hi = U +Uext( r i))- V ee stands for 
the electron-electron repulsion as a local operator. The 
exchange component of the perturbed operator is taken 
into account by 

occ occ 

v OEP (r) = ^if» = - k(r)}, (27) 



where Vh(t) = ^2° cc ji(r) and fcj(r) are the correspond- 
ing local Coulomb and exchange operators, fcj(r) is given 
according to Eq. (15), where v^ EP (r) = Ym=i ^( r )- The 
partition of the Hamiltonian given above can be taken 
as an operator for M0ller-Plesset (MP) perturbation ex- 
pansion [^9],[30| employed on orbitals and obtained by the 
solution of x-only OEP equations. 

The exchange-only OEP total energy Ef® EP can be 
considered as the sum of zeroth and first-order energies, 



E. 



•xOEP 
tot 



E^+E^=J2^ -<*« |W|* fl ) (28) 

i 

occ 



where 



^ occ „ 

E^ p [{u i }] = -\Y,jdTdv 



, Ui (r)u*(rK(r'K-(r') 



(29) 



According to Eq. (9) the exact (correlation corrected) 
total energy E^' OEP can be given at the nth order level 
of perturbation theory 

4 " ) ' OBP [{-J] = ^ ? £P + ^ 2) . (30) 

The first correction to the exchange-only energy occurs 
in the second order of perturbation theory 
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The correlation energy is expanded as the sum of cor- 
relation contributions at any order, so that 



4 n) [{ui>ei}] 



p=2 



'[{^i.Ci}]. 



(31) 



Next step is to express the unknown correlation po- 
tential in order to solve the correlated KS problem given 
by Eq. ( |22] ) . The correlation potential is also determined 
according to the perturbation theory at nth order level 
of theory keeping in mind that correlation and exchange 
effects are decomposed, 



v xc (r)=v x (r)+v^(r). 



(32) 



For v c we employ the expression, obtained from the exact 
transformation of the OEP integral Eq. (jlj) and given 
by Eq. (15). For sake of simplicity, however, we use 
the remarkably accurate KLI approximation, and neglect 
the rather difficult last term in Eq. (15). According to 



Eq. m 



can then be written 



The nth-order orbital dependent potential can be 



given using the M0ller-Plesset correlation energy E, 
given in Eq. ( [jl| ) 



(n) 
V ■ 



(r) = 



(n) 



1 SE, 



u* 5uj 



(38) 



In the next section we will show that the functional 
derivation with respect to the orbitals can be given ana- 
litically and finally one can get closed form for the cor- 
relation potential at nth-order level of PT. 

According to Eqs. (22-26) the following self-consistent 
Kohn-Sham procedure can then be constructed, which is 
correlated at the nth order level of RS PT, 



1 



-V 2 + v ext (r) + v h (t) + v x (r) + »W(r)]uj(r) (39) 

= e»u»(r)- 

We will discuss the various properties of this self- 
consistent KS PT scheme in section V. 



,(«) 



(r) 



(r) 



ij 



p(r) 3 3 



v { :h (33) 



We would like to note here that none of the existing DFT- 
PT schemes J^P^ / used directly this OEP expression 
for the correlation potential. Instead they derived the 
correlation potential in an alternative way. Gorling and 
Levy employed functional derivation over the second or- 
der correlation energy with respect to the Kohn-Sham 
potential and to the eigenvalues. Engel et al. also ap- 
plied the GL PT over the relativistic OEP leading to 
functional derivations with respect to the KS orbitals and 
to the KS eigenvalues. Both theories will lead to different 
expressions. We use, however, functional derivation with 
respect to only the orbitals according to Eq. (13). In 
the next few steps and in the next section we will show 
that this difference will lead naturally to different correla- 
tion potential then those given by the authors mentioned 
above. This formulation of the correlation potential has 
Slater and response component ||||, 



,(«) 



(r) 



v. 



resp,(n) 



and can be extracted as follows, 



(34) 



IV. THE CORRELATION ENERGY AND 
POTENTIAL AT SECOND ORDER 

Because of the increasing complexity of perturbation 
expansion, we give here only the MP perturbation energy 
at second order. For a general case one can give (in prin- 
ciple) those quantities at higher order as well. Note that 
one has to consider not only double excitations but also 
single configurations. In case of HF theory £?( 2 )> SOT ^ e = q 
by Brillouin theorem, but in case of KS self-consistent or- 
bitals we still have to calculate it |30|-^2]| . This term does 
not vanish completely but must be rather small for KS 
orbitals as well. For a closed-shell system, the second- 
order energy can be written in terms of sums over spatial 
orbitals as 



Ef\{u i }) = E l 



2, single 



E. 



2,D 



E 



\($ s \V ee {r)-V U ^(r)\$ s , r )\ 



EE 



r=l 

\w ik \ 2 

Ej - £fe 



E„ — E, 



EE 

ij kl 



| W «|2 



e k - ei 



(40) 



;fW(r) = 



OCC f s 

p(r) C3 

3=1 PV ' 



'(r), 



(35) 



j 



v 



S(n) = 



drp 3 -(r)«fW(r), 



(36) 



-OEP 



(i\v2 EP (r)\k), 



-M 

v ■ 



drpj(r)v^' (i 



(37) 



w^H$ s |v ee |^) = |<u||fcOI 2 , 



(41) 
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where |(zj||fc^)| 2 

(ij\kl) = I dvdt 



(ij\kl)(2(kl\ij) - (kl\ji)), 
,Ui(r)uj(r')uk(r)ui(r') 



v?f>W([p,{u i ,e i }]-,r)=v s x (r)+v s c > D >W(r) 



(47) 



r — r 



(42) 



occ vir 



mi? 



ijkl 



ij M 



$ s is the TV-electron Kohn-Sham ground-state single 
determinent wave function. The rth excited state to 
the Kohn-Sham equation is given by $ s . r . The ex- 
pectation values are as follows, E s = (& S \H \& S ) and 
E s .r — (®s,r\H \& s ,r)- Note that since the quantity be- 
ing summed in Eq. (^) is symmetric in i,j,k and I it 
vanishes when i = j or k = I. 

The next question to be addressed is the derivation of 
the correlation potential v c according to Eqs. (15) and 
(18) bearing in mind the decomposition v xc — v x + v c . 
The correlation potential at a given order of the pertur- 
bation series v c n ^ (r) can be given formally by means of a 

functional derivative of the correlation energy E c with 
respect to the orbitals Ui(r) according to Eqs. Jl3| ) and 
(38) via the orbital dependent quantity 



and the generalized response part 
< e / pA(2) ([P,{^}];r) = v r x esp (r)+v r c e ^ (2) (r) - 

m—1 s occ „ occ 1 

E 7 A ^ { E/ ^Efo-N - \{ij\U)) 

mi? 



ikl 
occ vir 



i kl 



EE / *tE 



mm? 



(i + £j — (k — (I 



(48) 



1 



(2) 1 SE^[{ Ui }} 

Vci 1 ' u*(r) S Ui (r) 

g single ^(2),D 



(43) 



+ 



u*(r) { Sui(r) Sui(r) 



According to Eqs. (38) and (41) the quantity vj?' at 
second order 



(2) / \ ( > ..-w i: fill I ^ I .1/ , I i \ 

«ci( r ) = u ci + v d = ( 44 ) 

Svi 



(2), single (2),D 



EE 

i fc>i 



[?I«] S 



Ei + e, - e k - ei 



The expression is given only up to the second order, and 
the extension of the formula to include higher order terms 
is straightforward , although order-by-order with increas- 
ing complexity |3(J . 



V. THE SELF-CONSISTENT PERTURBATION 
SCHEME 

Having discussed the partition scheme and the formu- 
lation of nth-order xc-potential, we would like to study 
the procedure given by Eq. (39). This scheme allows 
one to improve the eigenfunctions and eigenvalues of the 
unperturbed problem systematically. The nth-order cor- 
relation potential v < c n \r) (Eq. (|33|)), which is calculated 
on top of the first-order level of perturbation theory, 
is substituted back to the exchange-only KS equation 

(n) 



mi? 



5ui(r) 



\(ij\\kl)\ 2 = \j\kl](2(kl\ij) - (kl\ji)) (45) 



-[kl\j](ij\kl) 



\j\kl] = Mr) / dt 



jUj {r')ui{v') 



The Slater component of the second order correlation 
potential given by Eq. (|35|) will be precisely 



„£. «(r) = v S^gle,(2) {r) + y S,D,{2) _ 



(46) 



Using the decomposition v xc = v xc + v r xc sp , the potential 
part of v xc or equivalently the generalized Slater com- 
ponent of Eq. (|l8|) at second order of PT can be finally 
written (note that only double excitations are included, 
^smgie,(2) ^ g gj yen j n Appendix), 



in order to achieve self-consistency in the field of v. 
The introduction of Vc via Eq. ( |33|) together with Eq. 
(39) opens a new variational freedom in the original KS 
Eq. (E), which leads to a new iterative scheme. A new 

(n) 

v c can then also be constructed on top of the variation- 
ally solved correlated KS-problem, however, this would 
lead to a more complicated iterative scheme. 

The iterative scheme given by Eq. (39) is working in 
the following way: 

(i) one has to solve the first order level of theory which 
is exactly the exchange-only OEP reference state in this 
particular case providing the set of orbitals and eigenval- 
ues, {ui} and {ei}. 

(ii) the next step is the construction of the nth- 
order correlation energy (Eq. (|3l"|)) and the exchange- 
correlation potential v x n ,} on top of the exchange-only 
reference state (Eqs. ( |l8| ) and (|33|)). 

(iii) substituting back the potential Vc obtained in 
step (ii) to the KS equation a new correlated KS problem 
can be solved. In this way a new set of eigensolution 
{ui,€i} can be obtained. One can stop at this point, 
when the nth-order correlated KS problem is solved, and 
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the ground state energy can then be calculated by Eq. 
(9). 

(iv) Although present paper is addressed to the formu- 
lation of a standard perturbation theory on top of x-only 
OEP (steps i-iii), we indicate that further iterative steps 
can also be considered in practice when n ^ oo. The 
orbitals and the potential v x obtained in step (hi) 
differ from those obtained in step (i). Therefore, a new 

E ( c n) and i4™' 1 can be calculated which, however, leads 
to a KS problem again. This procedure continues until 
convergency is reached, provided, of course, that such a 
method turns out to be convergent. These additional it- 
erative steps with the "updated" v c n ^ result, however, in 
the change of the unperturbed problem. This would lead 
to the change of the partition scheme, the formula given 
by Eq. (40) for E c n ^ e.g. at second order of PT. There- 
fore, in Eq. (40) v (r) contains not only v x but also 

v c 2 ^ (r) so that, 

M/( 2 ) ( r ) = V ee (r) - v OEP (r) + vf ] (r). (49) 

Here we restrict ourself to the discussion of iterative 
procedure with only unchanged partition scheme of the 
Hamiltonian. The more comprehensive discussion of such 
a RS scheme with improved partition of the Hamiltonian 
can be the subject of further studies. 

The solutions of the correlated KS equation obtained 
in step (hi) can only be taken as final solutions of Eq. (39) 

(n) 

when infinite order PT is employed for calculating v c 
(n = oo). This is because the RS theory is mathemat- 
icaly equivalent representation of the full configuration 
interaction (full CI) theory |30| at infinite order expan- 
sion. In practice, only truncated perturbation expansions 
can only be considered, and, as such, fast convergence in 
the PT expansion is of great importance. Therefore, the 
proper choice of the reference space is also important to 
achieve rapid convergence |Q . Note that in this scheme 
the reference state is subsequently optimized when steps 
(i)-(iv) are used. 

Like in standard Hartree-Fock or KS procedures, the 
single particle Eq. (39) have to be solved iteratively un- 
til self-consistency is reached. According to Eq. ( |33| ) the 
scheme can be improved systematically by considering 

(n) • .rr 

higher order v c . One of the main differences between 
the original form of the RS scheme and which is given in 
Eq. (39) is that in the OEP PT scheme the wave-function 
and the eigenvalues are improved via the solution of the 
nth-order KS problem. This theory never goes beyond 
the single-determinant picture in the present form. Ex- 
ctited configurations are considered, however, when cal- 
culating Ec and Vc ■ This scheme alows one to drive 
a system progressively from a noninteracting reference 
state toward a fully correlated system. When step (iv) 
is switched on, the orbitals {u.i} and eigenvalues {e^} in 
the RS expansion of E^ (Eq. (|3l|)) are allowed to be 
the updated ones in the RS perturbation expansion. The 
scheme can therefore be considered as a self-consistent 



RS perturbation problem (SC-RS OEP). An alternative 
approach to carry out PT calculation sclf-consistently is 
the variation-perturbation approach [jl8|,[54j. Within the 
exact SC-RS OEP scheme self-consistency can be reached 
once the exact v xc is evaluated. We hope that this itera- 
tive scheme might have the " charming feature" of includ- 
ing the important part of electron correlation already at 
the second order of the correlation energy. Normally, the 
eigenvalue problem up to the first order level of theory is 
solved variationally. The iterative scheme (39) permits 
getting the best optimized first order reference state in 
RS theory. The optimization of the reference state is ac- 
complished by the reiteration of the updated xc-potential 
and by the consequent solutions of the "updated" KS 
problems until convergency is reached in step (iv). The 
self-consistent scheme given by Eq. (39) is similar to the 
so-called density functional perturbation theory (DFPT) 
]l8[ in that respect that succesive orders of perturbation 
are obtained iteratively in both cases. However, DFPT 
uses energy derivatives with respect to the ordering pa- 
rameters as a variation-perturbation treatment. 

Gorling and Levy (GL) Jll|,^8| have given more gen- 
eral functional for E xc and for v xc on the basis of PT. 
They found on KS basis that both the xc-energy and the 
xc-potential must be complicated functionals of both or- 
bitals and eigenvalues and also of a linear response func- 
tion type inverse operator G _1 . However, the quantity 
G" 1 is analytically unknown. The GL perturbation the- 
ory can provide an exact formal Kohn-Sham scheme only 
in basis set representation. Both the GL and the OEP 
PT theories lead to the exact formal representation of 
the KS equation, however, the construction of the nth- 
order v xc is different. In the GL PT a coupling constant 
a dependent scheme links the noninteracting A-electron 
system with the interacting real system where the elec- 
tron density remains independent of a pi| . The correla- 
tion potential and the energy are expressed in a Talyor 
series with respect to the a. The OEP PT scheme is a 
coupling-constant free formalism and the electron density 
does not remain constant during the perturbation treat- 
ment. The exchange-only OEP will result in different 
density from Eq. (39), although the difference may be 
rather small. In OEP PT a particular form of v xc is used 
given by Eq. (15), while in GL PT a more general, but 
unknown form of v xc is used. The main advantage of the 
present scheme is that the functional derivative Eq. (|3^ ) 

is directly accesible and, therefore, the nth-order poten- 

(n) 

tial Vc can be given analitically. Both theories deliver 
the exact exchange-correlation energy and potential or- 
der by order. Holas and March give the leading term 
of GL PT correlation potential in terms of first and sec- 
ond order density matrixes (2DM). The 2DM is also ex- 
pressed by PT. They also derived an exchange-potential 
which is free of the energy denominators ti — €k while 
the corresponding expression in GL theory does contain 
it . The extension of relativistic exchange-only OEP 
with perturbation theory has first been suggested by En- 
gel et al. very recently Although their scheme is 
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closely related to our OEP-PT approach, the construc- 
tion of the correlation potential seems to be somewhat 
different. Not taking into account its relativistic features 
functional derivatives appear in their approach with re- 
spect to the eigenvalues which are not included in the 
present scheme. We use Eq. (33) for the correlation po- 
tential together with Eq. (38) which will lead to some- 
what differnet formula than that given by Engel et at. 
Although they give a very accurate and detailed formu- 
lation of the theory we belive that our approach is still 
relevant for applications in the next future. 

In our SC OEP PT scheme the exchange-correlation 
potential can be given according to OEP method. The 
most convenient way is to use the KLI approximation |jj 
for the xc-potential according to Eq. (|l8[). The exchange- 
only KLI approximation given by Eq. ( |l8[ ) performs re- 
markably well |3]]l4|| providing nice agreement with the 
Hartree-Fock total energies. In the GL scheme one has a 
relatively complicated formula at the first order level of 
PT, e.g. at the exchange-only level using orbital, eigen- 
value and G _1 dependent exchange-potential. The SC 
OEP-PT method uses the much simpler exchange poten- 
tial according to Eqs. (15) and ( p9|) which is only orbital- 
functional. However, the OEP PT scheme given in Eq. 
(39) leads to completely exact theory when the last term 
in Eq. (15), which is neglected in the KLI approxima- 
tion, is calculated as well. Although the integrodiffer- 
ential equation Eq. ( |l2[) can be solved only numerically, 
however, its exact transformation given by Eq. (15) can 
be computed exactly for arbitrary fermionic system. The 
neglected term in Eq. (15), when considered, can increase 
the computational difficulties significantly. 



E 



in) 



is also the functional of the optimized exchange- 

(„) 

XC 



correlation potential v- 
by Eq. (39), 



via the iterative scheme given 



4? = 4? 4?] 



(50) 



This is because in steps (hi) and (iv) in Eq. (39) the cor- 
related KS problem is solved using the "updated" Vc . 
This equation reflects the complexity of the eigenvalue 
problem given by Eq. (38). Also, Eq. (|5p| ) demonstrates 
the difference of the SC-RS OEP scheme from the origi- 
nal formulation of the RS perturbation theory. The ex- 
panded nth order correlation energy is explicit functional 
not only of the self-consistent orbitals and the eigenval- 
ues, but also of the self-consistent exchange-correlation 
potential. Similar results are reported on Sham-Schlutcr 
basis 13|l4|] and by Engel et al. pl|. In that and other 
studies the variational energy expression has been treated 
as the functional of the fully dressed Green function G 
and of the reference function Go , and they are associated 
with the Dyson equation jjjfl. 



VI. CONCLUSIONS 

In this paper we have introduced a first-principle, 
parameter-free perturbational density functional scheme 
in which all exchange-correlation effects are consistently 
represented in terms of the eigensolutions of the Kohn- 
Sham equations. A self-consistent formulation of pertur- 
bation theory is developed on top of the exchange-only 
optimized effective potential (OEP) method. This gener- 
alization of OEP opens a new variational freedom which 
leads to a new iterative procedure. The total energy is 
also the functional of the nth-order exchange-correlation 
potential. In this scheme, the correlation energy and the 
correlation potential is expressed via perturbation series 
while the orbitals and the eigenvalues are variationally 
optimized via the Kohn-Sham equation. We give exactly 
the potential and the response parts of the exchange- 
correlation potential which are explicitly functionals of 
the eigensolutions of the Kohn-Sham eigenvalue problem. 
The correlation potential is given directly as a functional 
derivative over the M0ller-Plesset correlation energy ex- 
pression with respect to the Kohn-Sham orbitals. 

We also discuss the difference and similarities between 
the present and other perturbation theories which are 
also based on Kohn-Sham orbitals. 

Further investigations will be fruitful starting from 
the formalism introduced concerning the selection of the 
most appropriate reference state, the Hamiltonian par- 
tition scheme and the convergence of the RS PT based 
on OEP. We would like to emphasize the utility of Eq. 
(39) as a way of improving Kohn-Sham eigenvalues for 
calculations of band gaps, ionization potentials and ex- 
citation energies. The self-consistent OEP PT scheme 
introduced here can be useful in a wide range of areas in 
quantum chemistry and in solid state physics. Whether 
the present approach can be really competitive among 
known high accuracy but costly approaches, such as con- 
figuration interaction (CI), will be hopefully the subject 
of further tests in the future. 
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APPENDIX 

In order to get the second-order correlation potential 
one has to derive the single excitation components which 
occur when we are dealing with a " one-body" single par- 
ticle Hamiltonian. In Hartree-Fock theory single particle 
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excitations do not play a role in the second order PT cor- 
relation energy formula because (<I>, i |IT|<E , g fc ) = where 
$l k is the single excited determinent wave function. How- 
ever, its magnitude and contribution to the second-order 
correlation energy in PT-DFT must be rather small as 
well. Further arguments in this respect can be given only 
computationally. 

(2) 

Let us again decompose v c into contributions from 
single and double excitations, 

vi 2 \r)=vr^(r)+v^ 2 Hr). (51) 

v^'^ has been given already in section IV in Eqs. (47) 

and (48). Together with Eqs. (33)- (38) v ° ingle 'W can be 
partitioned to Slater and response parts 

v single,(2) _ y S, single, (2) _|_ y resp, single, (2) 



V 



(2) , single 



1 SE, 



(2) , single 



where 



g£,(2), single 

dm 



OEP , 



5v?ff SvLu 6v 



J x,ik 



■ — resp 
x,ik , u ^x,ik 



Sui 



(52) 
(53) 



ti-v 



mr z° cc (m} + {v\ k })^ — 



■OEP 
k 



(54) 
(55) 

(56) 



Sv r x esp 



m— 1 r 



P ^ 

F 3=1 



j, m— 1 



SAY,, 1 








UU X] 






5ui 



■HF\ 



(62) 



Sv 



xj 



3=1 



HF 



Sui 



(63) 



The functional derivative over the HF exchange energy 
density 



5e 



HF 



occ „ 
3=1 J 



,u*( r> 3 (i>*(r') 
|r-r'| 



(64) 



For the inverse matrix the following equality holds, 

(65) 



where 



Sui 



8u, 



3=1 



iPi 



(66) 



For simplicity, we give here the terms in conjuction with 
the remarkable KLI approximation || , therefore the last 
term given in Eg. (15) is neglected. However, this term 
is rather small ||[l4| and can be accounted for when high 
numerical accuracy is required. 



= «f (r)« fc (r) + Ui (r)u k (r) 6 -^-, (57) 



Sv x 

Sui p 



* a HF 



1 / u*e x 



Se 



HF 



Sui 



(58) 



where formally the HF exchange energy density is the 
one given in Eq. (E9) e.g., 



E 



OEP 



dre» F [{u° EP (r)}}, 
= 4 (r)p(r). 



6u 



J) V res P 

" A --v x ^{v)u k (r)+ Ul (v)u k {v)- ■■' 



Su 



(59) 
(60) 
(61) 
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